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Abstract 

The structure and the dynamics of homogeneous turbulence are modified by the presence of body forces such that the CorioUs or the 
buoyancy forces, which may render a wide range of turbulence scales anisotropic. The corresponding statistical characterization of 
such efl'ects is done in physical space using structure functions, as well as in spectral space with spectra of two-point correlations, 
providing two complementary viewpoints. In this framework, second-order and third-order structure functions are put in parallel 
with spectra of two-point second- and third-order velocity correlation functions, using passage relations. Such relations apply 
in the isotropic case, or for isotropically averaged statistics, which, however, do not reflect the actual more complex structure of 
anisotropic turbulence submitted to rotation or stratification. This complexity is demonstrated in this paper by orientation-dependent 
energy and energy transfer spectra produced in both cases by means of a two-point statistical model for axisynnmetric turbulence. 

^ We show that, to date, the anisotropic formalism used in the spectral transfer statistics is especially well-suited to analyze the refined 
5h dynamics of anisotropic homogeneous turbulence, and that it can help in the analysis of isotropicaUy computed third-order structure 

Q function statistics often used to characterize anisotropic contexts. 

^ Keywords: Homogeneous turbulence. Anisotropy. Spectra. Structure functions. 
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1. Introduction 

Homogeneous turbulence submitted to distortions such as 
solid body rotation, stratification, or the Lorentz force in the 
MHD context, exhibit axisynnmetric statistics, as a clear depar- 
ture from isotropy. These anisotropic efi'ects that arise due to 
modified dynamics or energy exchange do not fall within the 
classical description of turbulence in Kolmogorov's theory. The 
foUowing general questions, pertaining to the understanding of 
anisotropic homogeneous turbulence, can therefore be raised: 

1. How can we characterize the anisotropy of the flow, at 
one or two point, in physical or spectral statistical descrip- 
tions? 

2. Are passage relations available between: (a) the two-point 
statistics of structure functions, necessarily measured in 
physical space in experiments, and (b) kinetic energy 



and transfer spectra? (Which can be obtained for high 
Reynolds number turbulence from statistical models.) 

3. How can we interpret the modified dynamics of 
anisotropic turbulence? (It can be observed in terms of 
anomalous scalings of energy spectra or energy transfer 
spectra readily available in two-point statistical models.) 

4. Can we exploit some features easily accessed in spectral 
space to refine the phenomenological analysis of the dy- 
namics of turbulence in physical space? 

5. Do Kolmogorov scalings available for high Reynolds 
number isotropic turbulence apply to isotropically inte- 
grated statistics of anisotropic turbulence? 

6. What level of complexity do we have to introduce to go 
beyond the mere isotropic description, e.g. starting with 
axisymmetric flows? 
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As a bootstrap, we start, hereafter in this introductory sec- 
tion, by reviewing a few of the existing works related to these 
issues, and introduce some of the existing resuhs and formal- 
ism available in the statistical characterization of axisymmetric 
turbulence. In section |2] we describe the relationships avail- 
able for isotropic turbulence for second- and third-order statis- 
tics appearing in both Lin's and Karman-Howarth's equations, 
and suggest to extend the description of the velocity increment 
statistics to the anisotropic case, and to relate them to the mod- 
ified dynamics in axisymmetric homogeneous turbulence, es- 
pecially energy and transfer spectra. Results for the latter, in 
the case of stably stratified or rotating turbulence, are provided 
in section 4. 1 We choose to use an anisotropic two-point sta- 



tistical model, described shortly in section [3] since it permits 
to reach higher Reynolds number and smoother statistics than 
direct numerical simulation. Assuming isotropized passage re- 
lations, we then derive second- and third-order velocity struc- 
ture functions, presented in section [42} and discuss the results 
against existing scalings for isotropic turbulence. 

The spectral formalism for homogeneous turbulence allows 
to remove altogether the role of pressure, and hence permits a 
highly refined analysis of the dynamics of anisotropic turbu- 
lence, and on its sources; for instance from the role of inertial 
waves interactions in rotating turbulence, or for the dynamics 
of stably stratified turbulence in which the transfers are split be- 
tween internal waves and potential vorticity interactions. Even 
in flows which are thought to be isotropic, in experimental or 
numerical realizations, a degree of anisotropy may be hidden, 
depending on how the characterization of the isotropy of the 
flow has been done: the isotropy of nns velocity components 
or of integral length scales only characterizes the large scales, 
whereas vorticity or dissipation can be used for the smaller 
ones. The computation of anisotropic spectra allows to quan- 
tify such scale-dependent level of isotropy, or anisotropy, and 
to isolate the most important energy transfer contributions. Al- 
ternately, the statistics in physical space do not easily permit 
the removal of pressure from balance equations, and thus pre- 
vent from a tractable access to modal decompositions (e.g. the 



toroidal/poloidal decomposition). It however applies to locally 
inhomogeneous flows, which is not easily accessible to spec- 
tral analysis. We therefore emphasize the importance of both 
spectral and physical viewpoints, and the fact that they are com- 
plementary, which we expose and discuss in this paper, taking 
rotating or stably stratified turbulence as supporting evidences. 

In what follows in this section, we review a few studies which 
have been devoted to the characterization of anisotropy in tur- 
bulence, from the point of view of directional statistics in phys- 
ical space — including both anisotropy of velocity components 
or dependence on the separation vector — and studies about the 
anisotropic spectral scalings in anisotropic turbulence, related 
to the dynamics of the flow. 

If initially isotropic, homogeneous turbulence can be ren- 
dered anisotropic by introducing an external distortion on the 
flow: solid body rotation is present for instance in geophysical 
flows and acts through the Coriolis force, as well as a buoyancy 
force in density- or temperature- stratified layers. In conduct- 
ing fluids, the action of an external magnetic field also modifies 
the symmetries by means of the Lorentz force. In these three 
examples, the intensity of the corresponding force depends on 
the orientation of the fluid motion with respect to either the 
rotation axis, the gravity axis, the background magnetic field 
axis. Wave propagation can also be present, e.g. inertio-gravity 
waves or Alfven waves. They provide an anisotropic way of 
redistributing energy in terms of scale and direction, such that 
the dynamics of energy exchange in the turbulence is strongly 
modified. We focus here on the eff'ects of solid body rotation 
and stable density stratification on homogeneous turbulence. 

Studies about anisotropic turbulence were proposed and 
stemmed from various preoccupations. The theory of axisym- 
metric turbulence was regarded as a logical extension to the the- 
ory of isotropic turbulence, and was developed by Batchelor [I] 
— in which parallel and perpendicular directions (with respect 
to the axis of symmetry) were distinguished to express statisti- 
cal quantities, e.g. the dissipation — and by Chandrasekhar Q 
in the 50s. In these studies, the imposed symmetries do not ac- 
count for the possible presence of helicity, unlike the extension 



proposed by Lindborg about the kinematics of axisymmetric 
turbulence f3l. More recently, some exact vectorial laws were 
also proposed by Galtier for rotating homogeneous turbulence, 
exposing the need for a transverse/longitudinal components dis- 
tinction (with respect to the separation vector) when comput- 
ing velocity structure functions Q. An extensive discussion 
and review of the anisotropy in turbulent flows was proposed 
by Biferale & Procaccia Q using the symmetry group analy- 
sis (SO(3) for Navier-Stokes equations). These authors discuss 
Kolmogorov's theory and how to relate anisotropic flows to the 
«-th order structure functions, especially looking at the 4/5 law 
for the third-order one, which includes the isotropy hypothesis. 

On the one hand, the structure-function approach has re- 
cently been used for an experimental characterization of ro- 
tating turbulence dynamics by Lamriben et al. |6|, and for 
mesoscale turbulence in the atmosphere by Lindborg & Cho 
ITTI . both works including a discussion of the level of anisotropy 
to account for, as a departure to classical Kolmogorov scaUngs 
of velocity increments moments. In an experimental study con- 
cerning a turbulent jet, Xu & Antonia ||8J emphasize the im- 
portance of discriminating between the longitudinal and trans- 
verse velocity components when expressing the structure func- 
tion in axisymmetric turbulence, and Oyewola et al. use struc- 
ture functions to describe the anisotropy of the small scales in a 
turbulent boundary layer |(9J. In situ experiments by Kurien et 
al. in the atmospheric boundary layer also permitted to exhibit 
the different scalings of the structure functions, distinguishing 
between longitudinal and transverse increments, by separating 
the lowest order anisotropy contributions thanks to the SO(3) 
symmetry group. ifTOl 

On the other hand, several efforts were devoted to establish- 
ing the scalings that should apply to two-point velocity corre- 
lation spectra depending on the perpendicular ki_ or the par- 
allel wavenumber components with respect to the axis of 
symmetry. In the geophysical context, one may retrieve atmo- 
spheric spectra as in rotating strongly stratified turbulence. At 
small scales (several kilometers), Kolmogorov k']^^^ scaling is 
retrieved, whereas at large scales (several hundreds of kilome- 



ters), depending on the velocity component considered the scal- 
ing is either k"^ (zonal wind) or k'^ (meridional wind). IfTTI 
In the context of conducting fluids, plasmas and astrophysics, 
Galtier has also proposed several spectral scalings with spec- 
tra of the form k'^^k'^^^ to account for anisotropy in magnetohy- 
drodynamic turbulence fl2j. However, the physical arguments 
available for discussing the anisotropy of conducting fluid sub- 
mitted to the Lorentz force are not available for disentangling 
the intricate nonlinear dynamics of rotating turbulence, and, to 
some extent of stably stratified turbulence. 

In most of the above-mentioned works, spectra are used 
to characterize the cascade of energy and the dynamics of 
anisotropic turbulence, or structure functions compared with 
the isotropic theory scalings. Very few works are devoted to the 
relationship between the two formalisms, and the current paper 
is a first attempt at providing quantitative information on how 
anisotropic spectra dynamics can be used to interpret physical 
space velocity increment statistics. 

Thus, taking into account the modifications of the turbulence 
structure and dynamics due to external effects can be done at 
different levels. The statistical description of the velocity field 
can be done at a two-point level in physical space by the second- 
order velocity correlation tensor /?,j(r) - {ui{x)uj(x+r) where r 
is the separation vector Fourier-transforming this tensor yields 
the spectrum tensor 'I>ij(k) - 1/(8^)^ J Rij(r) exp(-ik ■ r)dr, 
whose trace is 0„(^:) = E(k)l2nk^ in isotropic turbulence, with 
k - \k\ and E{k) is the kinetic energy spectrum. Kolmogorov 
theory supports the scaling E{k) ~ k~^^^ in isotropic turbulence 
at high Reynolds number. 

However, external distortions render the turbulent flow statis- 
tics dependent on the orientation of r in physical space, e.g. for 
the two-point correlations Rijir). In isotropic turbulence, this 
tensor depends only on two scalar function, namely the longitu- 
dinal two-point correlation function /(r), and the transverse one 
g{r) fT3l, and is expanded over the tensors 6ij and rirj/r^. In 
axisymmetric turbulence, which is the case we consider here — 
where the axis of symmetry n is borne by the rotation vector 
in rotating turbulence or the axis of gravity in stably stratified 



turbulence — , Rijir) depends on four scalars and requires the 
addition of the two expansion tensors riiij/r and n,«j |1 1. 

The anisotropic dependence of the two-point correlation ten- 
sor in turn translates into a dependence of the corresponding 
spectra 'S>ij(k) on the orientation of the wavevector k with re- 
spect to the axis of symmetry n (the latter is the same as in 
physical space). The spectra can similarly also be decomposed 
in a general way over four scalars (see e.g. fT4\ and |fT5|). 

This correspondence between the physical space two-point 
statistical formalism of turbulence and the spectral description 
not only applies to second-order correlation tensors and spec- 
tra, but it also extends to third-order statistics, and moreover 
to velocity structure functions. The latter are of great use in 
Kolmogorov's description of isotropic turbulence, and are de- 
scribed e.g. in ||T6| and in IfTTll . Considering the velocity in- 
crement du - u(x + r) - u{x), in the isotropic context one 
commonly computes the «'''-order structure function {(6u)"} 
based on the longitudinal projection u of the velocity vector 
onto the separation vector r. The brackets indicate averaging 
using statistical homogeneity. Obviously, when turbulence is 
anisotropic, one should also distinguish different components in 
the structure functions and different orientations of r, although 
most of the theoretical results are available only in the isotropic 
context. One of them is the exact Kolmogorov four-fifth law 
wherein the third-order structure function Dmir) - {{6uf) - 
-4/5sr, valid for r in the inertial subrange, where s is the aver- 
aged dissipation ifTSl . or, in a different form proposed by Anto- 
nia et al. |fT9l : 



{6u(6uj6uj)} - --sr. 
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(1) 



where repeated indices imply summation. 

In the following, we investigate how anisotropic turbulence 
can recover the above scalings, using the two-point statistical 
EDQNM model. The question about the convergence of the 
statistics of the structure functions to the 4/5 limit has been ad- 
dressed before. Antonia & Burattini studied this limit at in- 
creasing Reynolds numbers Il20i . and it appears that the conver- 
gence is very slow, and may not be yet reached in the existing 



experiments. However, this low-Reynolds number effects has to 
be separated from the anisotropic effect in non isotropic flows. 
The leading order small-scales K41 scalings may be recovered 
in anisotropic flows, provided suitable filtering has been applied 
to the data, as shown with DNS fields by Taylor et al. 11211 . In 
the following, the interpretation of the scaling laws provided 
by EDQNM results therefore mixes the isotropic small-scales 
turbulent behavior with the anisotropic contributions, since no 
such separation as in Taylor et al. has been made. 

2. The link between spectral dynamics and velocity struc- 
ture functions 

The above-mentioned structure functions provide a scale- 
dependent statistical characterization of turbulence with a two- 
point separation. It is clear that they are related to the two- 
point correlation functions, and their spectral counterparts, 
such that second-order velocity structure functions are related 
to two-point velocity correlation spectra, and third-order struc- 
ture functions to energy transfer spectra (triple velocity correla- 
tions at two points). 

In the anisotropic context, and when dealing with statistically 
axisymmetric turbulence, the orientation of the Fourier vec- 
tor k with respect to the axis of symmetry borne by n, has to 
be taken into account, so that an additional dependence has to 
be introduced in the energy spectrum: E{k, ff). Integrating the 
latter over 9 e [-tt, tt] of course provides an equivalent spheri- 
cally averaged E{k). At this stage, the question of the resulting 
inertial scaling of the spectrum is raised. 

The kinetic energy spectrum E{k) and the nonlinear energy 
transfer spectrum T(k) appear in the Lin equation as 



d,E{k) + 2vt-E{k) = T{k) 



(2) 



in the isotropic case, where v is the viscosity, such that the to- 
tal dissipation is e = 2v J k^E(k)dk. In the rotating case for 
instance, phase-scrambling anisotropically changes the shape 
of the energy exchange term T(k, 6) in an equation equivalent 
to (j2]) (see next equation 



111. The transfer contains third-order 



correlation terms, and is thus related to the third-order statistics 



that are used to quantify nonlinearity such as the skewness of 
velocity gradient. 

Regarding second-order moments, one can obtain second- 
order structure functions starting with the two-point correlation 
spectra. For instance. 



functions 



Sin kr \ 

{{6ui6ui)) = 21 11 — E{k)Ak . 



Jo 



kr 
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where repeated indices are summed, contains information from 
all the velocity components. 

Equation (j2]) therefore relates the evolution of second- 
order statistics — the energy spectrum — to third-order ones — 
the transfer. The analogous in physical space of Q is the 
Karman-Howarth equation 

d, (m'V) = [dr + ^ j [RLL,L(r, t) + 2vdr {u'^f)] (4) 

in which f{r) is the longitudinal two-point correlation function; 
the longitudinal two-point third-order correlation function is 
Rhudr, t) = {uiix, t)u„t(x, t)ui{x+r, t)}r,„/r, r = |r|, and (3/2)u'^ 
is the total kinetic energy. Upon examining (|2]) and Q together, 
it is clear that if the dynamics of turbulence is anisotropic, one 
might expect not only a different equilibrium between the dissi- 
pative and nonlinear related terms of these equations, and a di- 
rectional dependence with k of the Lin equation for E{k) trans- 
posed into anisotropy with respect to r in the Karman-Howarth 
equation. The derivation of the latter equation relies on assump- 
tion of isotropy, but the general Karman-Howarth-Monin equa- 
tion, which contains the same terms that depend on the sepa- 
ration vector r can be derived as (see Frisch fVT], p. 78 for a 
derivation, and e.g. ||6|): 

l^{u(x)u(x + r) = ^V<5u(r)|5u(r)|2>+yV2<u(x)u(x + r>(5) 
2 ot 4 

The dependence of the velocity, or velocity increments 6u, 
on the direction of r with respect to n, is therefore an im- 
portant parameter in axisymmetric turbulence statistics. For 
isotropic turbulence, Kolmogorov's theory (1941) only con- 
siders a scalar longitudinal increment 6u, provides the scaling 
(((5m)") ~ (er)"^^, and allows to draw from equation ^ a sim- 
plified relationship between second- and third-order structure 



{{6uf} ^--^sr + 6vdr{{5uf). 



which yields, at infinite Reynolds number, the famous four 
fifths law 



= -(4/5)er, 



(6) 



or equation ([T]i. ifTOl 

A useful relationship can be estabUshed between third-order 
spectral statistics i.e. the nonlinear energy transfer T(k), and the 
third-order structure function : 



f 
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(6u {dqf) = 4r g(kr)T(k)dk 



(7) 



where g{kr) = (sin kr - kr cos kr)/(kr)^ and (6q)^ = (duidui). It 
shows the direct link of the third-order structure function on the 
dynamics of turbulence. 

Formulas similar to equations ([3]) and (|7]i rigorously derived 
for anisotropic turbulence are still not available in the general 
case, but, in the rotating case, the isotropic relationships ap- 
ply, since the Coriolis force is not a production term in the bal- 
ance equations, and drops out altogether in the isotropized bal- 
ance equations, both in physical and spectral spaces. In general 
anisotropic turbulence, the analogous of (|7]i implies a vector 
dependence of the structure functions on r. It also shows that 
a dynamics modified by an external distortion applied to ho- 
mogeneous turbulence translates immediately in a modification 
of the structure functions. It is therefore interesting to discuss 
the applicability of Kolmogorov scalings in flows that contain 
some anisotropy. Attempts at such scalings for rotating turbu- 
lence are proposed by Galtier ||4|, or by Lindborg & Cho using 
atmospheric data. [7| In the axisymmetric case, an extensive 
work is in progress Il22l . 

3. Statistical model for the spectral anisotropy of axisym- 
metric turbulence 

We study a simplified case which nonetheless captures some 
important anisotropic physical mechanisms in e.g. geophysical 
flows: stable stratification and rotation, taken into account in the 



Boussinesq system of equations for an incompressible fluid: 



d,u + u ■ Vu - vV H - -Vp - 2Q n X u + b n, (8) 



d,b + u-Vb-x'7^b = -N^nu, 



V H 



0, 



(9) 
(10) 



where b is the buoyancy field, associated with fluid density fluc- 
tuations around background density augmented by mean den- 
sity gradient. is the Brunt-Vaisala (buoyancy) frequency and 
D. the rotational frequency. The corresponding terms, Coriolis 
force and buoyancy force, act on the velocity field linearly and 
compete against the nonlinear advection term. In this study, the 
separate effects of either one is considered (although the combi- 
nation of both characterizes some atmospheric or oceanic flows, 
with a ratio a = 2Q./N -1/10). 

In rotating or stratified turbulence, the dynamics of the flow is 
superimposed with wave-like motion that transfers energy dif- 
ferently from the classical turbulent dynamics. In the linear 
limit, at large or Q, one recovers wave turbulence, as a soup 
of superimposed inertial waves or internal gravity waves. For 
instance, "turbulence" and "wave" like dynamics may be de- 
fined by splitting the velocity field in the eigenmodes of the lin- 



earized system ([8])-( 10 1, so that the total turbulent energy of a 
rotating and stratified flow can be divided into a "vortex" mode 
and a "wave" mode. The linear evolution of the wave mode is 
governed by he dispersion relation o-{k) = NsinO for internal 
waves and cr{k) = 2Q. cos for inertial waves, which depend on 
6, the polar angle with the vertical, whence a directional depen- 
dence of the turbulent motion. 

In order to study relatively high Reynolds number turbu- 
lence, we introduce a statistical model which describes the evo- 
lution of the two-point correlation spectra, whose dynamical 
equations derive from ([8|-([T0li, introducing a closure known as 
EDQNM. The Eddy-Damped Quasi-Normal Markovian model 
for isotropic turbulence was studied by several authors (see e.g. 
Orszag f23l) decades ago, although it has recently regained in- 
terest for its ability to reach very high Reynolds numbers as 
demonstrated by the spectrum we computed with EDQNM, 
shown on figure [T| (Similar figures are presented in ll24l .) 



The more advanced version developed for rotating or stably 
stratified turbulence (denoted EDQNM2 since it differs signif- 
icantly from the model for isotropic turbulence) reflects more 
accurately the wave dynamics and is capable of taking into ac- 
count anisotropic features. In the limit of very strong rotation, 
for instance, it becomes identical to the wave turbulence clo- 
sures described by Zakharov ||25J . although in this book, Za- 
kharov et al. do not consider anisotropic dispersion relations 
as for instance inertial waves. An extensive description of the 
EDQNM2 closure model is provided in 1261 fT4l . as well as com- 
parisons with Direct Numerical Simulations (DNS) and Large 
Eddy Simulations (LES) which show a very good agreement 
and illustrate the capability of the model to accurately repre- 
sent anisotropic turbulence. We find that, among other clo- 
sures, only the EDQNM2 model for closing the spectral trans- 
fer terms (the equivalent to T in equation |2]) was capable of 
predicting all the anisotropic features observed in the DNS and 
LES results. The model equations for the generalized transfer 
terms, e.g. T''(k) for the directional kinetic energy spectrum, 
involve sums of eight contributions (according to polarities of 
triadic interactions), weighted by the rotation-dependent factor 



(11) 



(or Brunt-Vaisala frequency in the stratified case): 

The numerator of the integrand takes into account the quasi- 
normal expansion for non-isotropic turbulence and is closed in 
terms of quadratic combinations of the two-point correlation 
spectra, whereas the denominator involves viscous and eddy- 
damping effects through 

^Tp/' = ^kpgil + 2\&kpgQ.(Ek\\/k + e'p\\/p + e"qn/q)) (12) 

for a triad of wavevectors k -1- p -H q - with polarities 
{e, e',e") = {±1)^, thus accounting for the explicit linear ro- 
tation effects on triple correlation through the phases. Thus, 
the classical timescale d-^pg in the isotropic non-rotating case is 
replaced by the preceding triadic complex timescales. As al- 
ready mentioned, the above expression is consistent with wave 
turbulence results |l22l|28l. The isotropic EDQNM model con- 
sists only of the Lin equations (j2]) in which the transfer term is 



closed using double products of energy spectra E(k), and where 
the only dependence variable is the wavenumber k. 

4. Second- and third-order structure functions derived 
from spectral statistics 

The EDQNM2 model for rotation and stable stratification 
consists of equations for the two-point correlation spectra sim- 
ilar to (j2]) in which the nonlinear transfer spectra in the right- 
hand-side are replaced by their closure. In the model for ro- 
tating turbulence, four coupled equations are solved for the en- 
ergy density spectrum e(k), the helicity /z(k) and the polariza- 
tion spectra Z(k), each of which is expressed at each discretized 
spectral point in the polar-spherical representation of wavevec- 
tors k, that is as functions of discretized (k, 6): 
d 



dt 



|-+2y/t2+4iQ^ 
at k 



+ 2vk^ e = 



JZ 



(13) 



dt 



+ 2vk^h = t'' 



The transfer spectra are also discretized in the same way. 

For stratified turbulence, the density spectrum of rescaled po- 
tential energy 03(k) and the cross-correlation spectrum for den- 
sity and velocity *P(k) are involved 



d_ 
dt 



2vr O 



^1 



kinetic/potential energy exchange, mediated by the flux *F pro- 
duced by buoyancy. 

The complexity of these systems of equations and anisotropic 
damping ( [T2| ) provided by the EDQNM2 model needs to be 
contrasted with the simplicity of the EDQNM model for E{k) 
whose dynamical equation is only (|2]l. 

The corresponding set of time-dependent differential equa- 
tions are advanced in time starting from initial conditions which 
correspond to a narrow-band distribution of kinetic energy cen- 
tered about a prescribed peak. The spectra then evolve into de- 
veloping an inertial subrange and a dissipative one. From this 
time on, turbulence evolves in a self-similar decay. 

The evolved rotating turbulence spectrum presented on fig- 
ure [T| corresponds to a Reynolds number Re = 7000 and the 
Rossby number is Ro = m/(QL) - 0.04 (u is the rms velocity, L 
the integral scale and Q. the rotation rate) although the figures at 
initial time were higher, since they usually decay in time. For 
the stratified turbulence spectrum. Re = 1000 and the Froude 
number is Fr - ujiNL) - 0.09, where is the Brunt- Vaisala 
frequency which characterizes the intensity of the vertical den- 
sity gradient. Note that the Reynolds number is lower in the 
stratified simulation, since a part of the dynamics is contained 
in the potential energy stored by density fluctuations around the 
mean density. 



4.1. Quick comments on the spectral statistics and the dynam- 



- +2vfe2|((]>2-H03) 



+ 2vr - 2iA^ sin 0k {(^2 - $3 -H I'Pr) 



(14) 



ics of anisotropic turbulence 
The developed spectra for isotropic, rotating, and stably strat- 



along with the kinetic energy spectrum decomposed in poloidal 
and toroidal spectra <I)2(k) and <l>i(k). 

Equations ( [T3| ) and ^\A\ are therefore generalized versions of 
the Lin equation Q and are exact equations, unless a closure 
such as equation ( [TT| i is applied. Moreover, in both systems 



(13 1 and (14 1, the decomposition of the spectra is optimized 
to yield simple and meaningful equations: the Z equation in 



( 1 3 1 contains the oscillatory part due to the effect of rotation. 



and the last equation of ( 14i represents oscillations due to the 



ified turbulence are plotted on figure [T| Note that spectra of 
anisotropic turbulence, in our model, but also in DNS or in 
experimental measurements, depend on the orientation of the 
wavevector, in agreement with two-point correlation functions 
which depend on the orientation of the separation vector The 
relevant parameter is the angle 9 with the axis of symmetry. 

We have purposefully chosen an very high Reynolds number 
for the isotropic spectrum for three reasons: (a) it demonstrates 
the capacity of statistical EDQNM closures to reach Reynolds 
numbers closer to those observed in geophysical or astrophys- 



ical contexts than is possible with DNS; (b) the anisotropic 
EDQNM2 closure requires additional computational effort, and 
therefore the inertial ranges observed on the corresponding 
spectra of figure [T] span roughly two decades and are reduced 
with respect to the isotropic case; however, equivalent DNS 
would require an order of magnitude larger computational ef- 
fort; (c) although the isotropic EDQNM spectrum exhibits six 
decades of inertial subrange, the corresponding kinetic energy 
transfer plotted on figure |2] only exhibits a narrow plateau at 
zero (if at all) between the energy-losing large scales and the 
energy-gaining small dissipative scales, also corresponding to 
constant downscale spectral energy flux. 

One should also note that statistical two-point models are 
by essence better adapted to produce smooth spectral statistics 
than Direct Numerical Simulations. First, the discretization of 
the spectral space can be adapted easily when solving EDQNM 
equations, by adjusting both the wavenumber k distribution and 
its orientations, with the possibility of accumulating grid points 
towards the equatorial or the polar directions. By construction, 
DNS is constrained to a Cartesian uniform grid discretization of 
Fourier components. Moreover, the recent experiment by Lam- 
riben et al. 16jj has shown that ensemble averaging requested 
to compute smooth third-order statistics — namely third-order 
structure functions or energy transfers — required hundreds of 
realizations, an issue which also applies to DNS realizations. 
Since the statistical two-point model equations already con- 
cern the statistics of the second- and third-order moments, ex- 
plicit ensemble averaging is self-contained in the model (as ob- 
served in the following smooth figures of results produced by 
EDQNM). 

The structure of decaying rotating turbulence has been de- 
scribed extensively from experimental ll29l . theoretical, and nu- 
merical results lfT4l[30l . Inertial wave turbulence was studied by 
Galtier ||311 . and Mirunni & Pouquet investigated energy and 
helicity spectra of rotating turbulence ||32]| . with diflterent spec- 
tral scaUng for each. The dynamics and modified energy cas- 
cade produce an anisotropic structure which was simulated with 
high resolution DNS also by Morinishi et al. Il33l . Anomalous 



scaling of structure functions in rotating turbulence also per- 
mitted Seiwert et al. to investigate the intermittency of rotating 
turbulence fM\. 

The interpretation of the dynamics leading to rotating homo- 
geneous turbulence structuration is still debated; the two main 
viewpoints rely either on linear timescales of inertial waves, 
with the Coriolis force acting on locally inhomogeneous struc- 
tures that emit waves |35|, or on a long-term nonlinear ef- 
fect fT4j. It is not the object of this paper to reconcile these 
viewpoints, but, in short, both agree that in rotating turbu- 
lence vortices are elongated along the axis of rotation. The 
corresponding dynamics is that of preferential energy trans- 
fer towards motion close to two-dimensional, although com- 
plete two-dimensionalization is not expected. It also means 
that energy accumulates in orientation dependent kinetic energy 
spectra E{k, 6) in the vicinity of the two-dimensional manifold 
9 ^ 7r/2 which is the neighborhood of horizontal wavevectors in 
spectral space. This is confirmed by the angular dependent en- 
ergy spectrum plotted on the top panel of figure [3] which also 
compares well with spectra processed from DNS of decaying 
rotating turbulence. Note that the power-law scaling of the cor- 
responding spectrum in figure [T] is the result of the spherical 
averaging over the directional spectra of figure |3] with no sin- 
gle identical power-law applying to each individual one. The 
bottom panel of figure [3] presents a tentative measure of the 
anisotropy in the spectra, by subtracting to the directional spec- 
tral the averaged spectrum. The curves present only the limit 
spectra, equatorial or polar; also note that the anisotropy is re- 
versed between the rotating and the stratified cases, with more 
energy in the polar direction (along the symmetry axis) in strat- 
ified turbulence, whereas the accumulation of spectral energy 
in rotating turbulence is in the vicinity of the equatorial (per- 
pendicular to the symmetry axis) direction. 

For the stably stratified case, the spectral structuration is re- 
versed, such that the energy accumulates in spectra of verti- 
cal k {6 =^ Q) as shown on figure [3] (see also e.g. Lindborg 
11 11 ). Again, the anisotropic statistical model compares ex- 
tremely well with experimental wind-tunnel observations and 



DNS ||26l. These comparisons not only concern dynamical 
quantities, such as energy or dissipation, but also directional 
integral length scales, and the most refined possible compari- 
son between directional spectra E(k, 0), that is for every direc- 
tion of every scale. In both the stratified and the rotating cases, 
EDQNM2 model predictions are very close to the same spectral 
statistics extracted from DNS fields, down to the most refined 
ones. In physical space, the structuration of stably stratified 
turbulence corresponds to a layering of the flow, with strong 
horizontal motion, and vanishing vertical velocity, although the 
flow retains a strong vertical variability due to the presence of 
large vertical gradients d/dz- 13611371 

The anisotropy observed in the spectra of figure |3] is the 
departure of each spectrum E(k, 9) around the average shown 
on figure [T] Since the plot is in logscale, the clear differ- 
ence between horizontal wavenumber spectra and vertical ones 
demonstrates the large ratio between the energy in vertical and 
horizontal motion. Moreover, one observes that in the strati- 
fied turbulence spectrum the anisotropy is largest at the larger, 
energy-containing, scales, although it extends to the middle of 
the inertial sub-range; in the rotating turbulence spectrum, the 
anisotropy extends throughout the inertial range down to the 
smallest scales. Of course, this observation depends on the 
values of the Rossby and Froude numbers as well as on the 
Reynolds number. The rotation (or stratification) timescale has 
to be compared to the local timescale of a given turbulent struc- 
ture associated with wavenumber k in order to assess whether 
its dynamics may be aff'ected by the Coriolis or the Boussinesq 
forces (for instance by computing a Zeman scale (Q^/e)^'^^ 
||38]| . where e is the dissipation, or similarly an Ozmidov scale 
for stratified turbulence by replacing Q by ||39l ). It nonethe- 
less demonstrates that there exist parameter ranges at which tur- 
bulence may be strongly affected by external distortions, and 
its structure and dynamics can depart significantly from that of 
isotropic turbulence. 

In addition to providing the second-order statistics corre- 
sponding to the two-point energy spectra discussed above, the 
closure provides quantitative information on third-order two- 
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Figure 1: Solid line: Kinetic energy spectrum E(k) of isotropic turbulence 
at Re^ = 26 X 10* given by tlie EDQNM closure model. Dot-dashed line: 
spherically accumulated kinetic energy spectrum given by the EDQNM2 model 
for rotating turbulence. Dashed line: spectrum given by the closure for stratified 
turbulence. The inertial range k^^^^ power-law is also shown on the figure, 
corresponding to the scaling of the isotropic spectrum. The spectra for rotating 
and stratified turbulence both scale like A-"' '. The wavenumber k is normalized 
by the inverse of the Kolmogorov length scale rj. 

point correlation spectra, i.e. the kinetic energy transfer spectra 
shown on figure|2] Apart from the Reynolds number difference, 
the spherically averaged transfer spectra of anisotropic turbu- 
lence are similar to the isotropic turbulence transfer. Note that 
we present spectra of the kinetic energy transfer T"{k), which 
is only a part of the energy transfers occurring in stratified tur- 
bulence, since there is also coupling with the potential energy 
mode. As observed in the isotropic case, in the anisotropic 
cases energy is drawn from the large scales and injected in 
the small scales in a classical forward cascade. In the mean 
time, if one observes the orientation-dependent transfer spectra 
(not presented here), one also notices a re-distribution of energy 
among different orientations of wavevectors, even at constant 
wavenumber, i.e. without interscale exchange. 

4.2. Structure functions 

The above observations of second- and third-order correla- 
tions spectra have a direct impact on the second- and third- 
order structure functions of turbulence, according to the dual- 
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Figure 2: Solid line: Kinetic energy transfer spectrum T{k) of isotropic tur- 
bulence given by the EDQNM closure model (scaled by 1 0^ to render it vis- 
ible). Dot-dashed line: spherically accumulated kinetic energy transfer spec- 
trum T'{k) given by the EDQNM2 model for rotating turbulence. Dashed line: 
transfer spectrum given by the closure for stratified turbulence. 

ity of spectral and physical space. Some passage relations be- 
tween the two spaces have been presented in section|2] We now 
present the results obtained by using these relations applied to 
the EDQNM closure results for isotropic and anisotropic turbu- 
lence. It is however not possible to derive any «'''-order struc- 
ture function from spectral information contained in two-point 
statistical closures, precisely because the spectral information is 
restricted to «'''-order correlations at two points. A case-by-case 
inspection is thus required. 

The relation (|3]l linking the kinetic energy spectrum to the 
second-order structure function is applied to the results ob- 
tained by EDQNM, and also to the anisotropic turbulence re- 
sults. The numerical resolution of the isotropic EDQNM model 
directly provides discretized distribution of E{k) and T{k). For 
the EDQNM2, the spherically averaged spectra E{k) and T{k) 
are obtained from the anisotropic spectral data E(k,6) and 
T''(k,6), by averaging over 0. Once these data are obtained, 
a simple numerical quadrature of the integrals in equations ^ 
and I?]) is used to obtain {(duidui)} and {6u {5qf) (likewise for 
{(Suf) with an integral not recalled here). The corresponding 
structure function is normalized by {erf'l^ and plotted against 
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Figure 3: Top panel: Orientation-dependent kinetic energy spectra E{k,6) ( 
seen to be monotonically varying with 6). Solid lines: results from EDQNM2 
model for stratified turbulence; dashes: results from the EDQNMt model for 
rotating turbulence, shifted down one decade. For each case, a subset of the 
discretized orientations 9 are presented. Spectra corresponding to horizontal 
wavevectors are indicated by black dots. Bottom panel: A measure of the 
anisotropy showing the difference between the polar and equatorial spectra with 
the spherically averaged one. Solid lines: EDQNM2 result for stratified turbu- 
lence, dashed lines for the rotating results. 
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rjrj {rj is the Kolmogorov length scale) on figure]?] showing 
the scaling which is obtained for isotropic turbulence, and il- 
lustrating the advantage brought by achieving a high value of 
the Reynolds number. For rotating turbulence and for stably 
stratified turbulence, the corresponding values of the Reynolds 
number are lower, and one has to deal with a narrower curve. It 
is clear, however, that the isotropic scaling is far from being ap- 
plicable to the anisotropic turbulence structure functions, since, 
here, if we compute max(((5M/5M,))/(£r)^^-'), one obtains 3.4 for 
the isotropic case, 4.2 for the stratified case, 6.6 for the rotating 
case. 

This is also observed on the curves for the third-order struc- 
ture functions plotted on figures]?] and ]6] in two different forms. 
Figure ]5] shows the averaged cubed longitudinal increment 
{(Su)^} divided by the behavior proposed in Kolmogorov the- 
ory: -{sr). The expected proportionality constant 4/5 is clearly 
recovered by the EDQNM closure transfer processed through 
equation (]6|, with a nice plateau tangential to the y - 4/5 hori- 
zontal line. It should be noted that even at this very high value 
of the Reynolds number the extent of the plateau is less than 
a decade in scale. This may explain part of the early contro- 
versies about the precise value of the proportionality constant 
obtained from experimental measurements. The same consid- 
erations apply for the other form of the third-order structure 
function proposed by |fT9l , {6u{6uj5uj)), with a proportional- 
ity constant 4/3 with er, plotted on figure]6]using equation (JTJ. 
In both plots, the scaling is obviously not the right one for ro- 
tating turbulence or for stratified turbulence, although the or- 
der of magnitude of the peaks of the structure functions for the 
anisotropic turbulence cases are consistent with the isotropic 
turbulence one. {•max{{6u{6uj6uj)))l{er)) is marginally lower 
than 0.8 in the isotropic case, 0.39 for stratified results and 1.02 
for the rotating results.) 

In rotating turbulence, the energy cascade is known to be 
reduced by the presence of rotation, so that the dissipation is 
lower than for an equivalent Reynolds number isotropic tur- 
bulence dynamics. This can explain the fact that the peaks 
of -{6u{6uj6uj)) lier) or of -{{Suf) lier) are higher than in 



isotropic turbulence. In addition, one observes on figures|5]and 
]6] a negative excursion of the third-order structure function in 
the very large scales, which could be attributed to a trend to- 
wards a reverse cascade, although this has to be confirmed by 
other simulations. 

In the stratified case, the trend is reversed, an eff'ect which 
could be attributed to augmented dissipation. The analysis is 
subtle here, since there is not only a cascade of energy due to 
classical nonUnear turbulent interaction, but there is also a cou- 
pUng between the kinetic energy and the potential energy. It 
seems that the balance between energy taken away from the ve- 
locity field for creating density fluctuations is more than com- 
pensated by the increased dissipation from the added scalar cas- 
cade. An attempt at including the part of dissipation arising in 
the potential energy balance fails at correcting enough the scal- 
ing of the third-order structure functions of stratified turbulence 
in figures]5]and]6] 

Finally, in order to address the Reynolds number dependence 
issue for the third-order structure function convergence, we 
computed an isotropic turbulence EDQNM result for a lower 
Reynolds number than proposed above, at Re - 2600, of the 
same order of magnitude as those for the anisotropic EDQNM2 
data. The solid curve of figure ]5] shows the corresponding re- 
sult, which underestimates significantly the 4/5 constant. The 
departure amplitude is similar to that of the anisotropic data, al- 
though we believe it is of different nature. Following the anal- 
ysis by Taylor et al. 11211 . there are therefore two sources of 
departure to the K41 scaling, one due to low Reynolds number 
effect, the other to anisotropy. Here, for the anisotropic runs, 
anisotropy is responsible for the departure from the asymptotic 
large Reynolds number scaling, since for rotating turbulence 
it leads to an overestimate, whereas for stratified turbulence it 
is underestimated. Also, the measures presented by Antonia & 
Burattini for isotropic turbulence exhibit a uniform convergence 
from below the value 4/5 EOl . Thus, although there may be a 
finite Reynolds number effect, we believe that, to first order, the 
difference with the high Reynolds number isotropic scaling is 
due to anisotropy. This argument is also supported by the curve 



for the second-order structure function in isotropic turbulence 
at the lower Reynolds number Re - 2600 presented on figure|4] 
It shows that although the asymptotic large Reynolds number 
scaling of is not yet reached, the maximum of {{SuiSu,)) j {eff-l^ 
is within 30% of the asymptotic value, whereas the peak for 
rotating turbulence is twice this value. The shape of the spec- 
tra and the difference in the downscale cascading rate in the 
anisotropic runs with respect to the isotropic one can only ex- 
plain this departure. 

5. Conclusion and perspectives 

We have described in this work the spectral statistics of 
anisotropic turbulence, with spectra of second- and third-order 
two-point correlations. These kinetic energy spectra and kinetic 
energy transfer spectra are the result of the spherical averaging 
of spectra of correlations taken with separations along differ- 
ent angles with respect to the axis of symmetry, or in other 
words for different orientations of the wave vectors. Not only 
these angle-dependent spectra correspond to variable intensities 
of the energy, but their inertial range scaling, when averaged, 
produces the power-law scaling of the overall kinetic energy 
spectrum. Results obtained by a two-point statistical closure of 
EDQNM type show that, at high Reynolds number, the result- 
ing anisotropy can be strong, and corresponds to a non clas- 
sical structuration of the flow. The cases of strongly rotating 
turbulence and strongly stratified turbulence were investigated, 
with respectively an elongation of the vortices along the rotation 
axis or a horizontal layering. The velocity components are con- 
sequently of variable magnitude in the vertical and horizontal 
directions, so that an impact on the velocity second-order struc- 
ture functions is expected. The access to the detailed energy 
transfers also permits the computation of third-order structure 
functions, and the comparison with exact scalings obtained for 
isotropic turbulence. As expected, the constants obtained are 
different, although thanks to the relationships presented in sec- 
tion|2]one is able to closely related the anisotropic spectral char- 
acterization to the structure functions, thus to link the dynamics 
of energy cascade and inter-orientation transfer to the statistical 



moments in physical space. The departure from K41 theory for 
anisotropic turbulence statistics of the second- and third-order 
structure functions shows that the Karman-Howarth equation 
from which theoretical scalings are derived exhibits a modified 
equilibrium when the structure of turbulence is anisotropic. The 
correspondence with the Lin equation shows that the dynam- 
ics of anisotropic turbulence is also modified, and we are able 
to trace the origin of the modified energy spectra scalings to 
the transfer terms. Since these transfer terms are accessible by 
spectral theory, this may be a way to improve the prediction of 
the third-order structure function moment in rotating or strati- 
fied turbulence, using two-point statistical closures. 

Of course, this calls for additional work and further devel- 
opments, including the derivation of relationships equivalent to 
([3]l and ([T} but considering the specific components of veloc- 
ity parallel or orthogonal to the axis of symmetry. In so do- 
ing, one will highlight the value of the anisotropic EDQNM2 
closure which not only allows to model high Reynolds num- 
ber turbulence but also provides quantitative information on the 
anisotropy of statistics. 
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